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Abstract. The pre-stellar cores in which low mass stars form are generally well magnetized. 
Our simulations show that early protostellar discs are massive and experience strong magnetic 
torques in the form of magnetic braking and protostellar outflows. Simulations of protostellar 
disk formation suggest that these torques are strong enough to suppress a rotationally supported 
structure from forming for near critical values of mass-to-flux. We demonstrate through the use 
of a 3D adaptive mesh refinement code - including cooling, sink particles and magnetic fields - 
that one produces transient 1000 AU discs while simultaneously generating large outflows which 
leave the core region, carrying away mass and angular momentum. Early inflow/outflow rates 
suggest that only a small fraction of the mass is lost in the initial magnetic tower/jet event. 
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1. Introduction 

The early evolution of low mass, isolated protostellar cores in the pre-Class stage is 
now understood as the interplay between gravity, rotation, magnetic fields, and radiative 
coolin g. Strong magnetic fluxes have been observed in molecular clouds (e.g. ICrutcher 
1999). Recent simulations have focused on the changes end ued by a magnetic field in 

the collapse, through either a 3D SPHR or AMR approach ( Ba.su fc Mouschoviasl 1994 ; 

Machida et al.lbOO^lHosking fc Whitwortll2004lMachida et al.l2005bUal;lBaneriee fc Pudritz 



20061 iPrice fc Batefl2007t iHennebelle fc Fromand 120081: iDuffin fc Pudritzj 120091) . It has 
been shown that magnetic fields slow the collapse timescale and brake the rotation of the 
initial core and of massive disc-like structures that subsequently form. Magnetic fields 
have also been shown to suppress fragmentation and the formation of bars and spiral 
waves. Furthermore, they facilitate the launching of molecular outflows. 

Several theoretical problems have arisen from these simulations. First, from axisym- 
metric simulations o utside of 6.7 AU, it is argued that Keplerian discs cannot form in an 
ideal MHD collapse (jMellon fc Lil 120081 ). except i n the limit of very weak m agnetic flux 
or high magnetic diffusivit y whether numeric al ( Krasnopolskv et al. 201fj() or through 
strong ambipolar diffusion ( Mellon fc Li 2009t ). Secondly, magnetic tension and pressure 
seem to efficiently suppress fragmentation in the early stages of protostellar collapse so 
that it is difficult to form multi-star systems. With the ability to extend the magnetized 
collapse further, we can begin to examine questions such a s how the Core Mass Fun ction 
(CMF) is related to the Initial Mass Function of stars fe.g. lMatzner fc McKeell2000l ). the 
nature of the molecular outflow on larger scales and whether fragmentation is suppressed 
even in later stages. We present our results of the early collapse using 3D ideal and 
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non-ideal magnetohydrodynamic simulations, and the first results of the later stages of 
the magnetized collapse and outflow. We are able to evolve the simulation an additional 
10 4 — 10 5 yr due to th e implementation of the sink particle in the FLASH AMR code 



( Federrath et al.ll2010l) 



2. Numerical Methods and Initial Conditions 

We model our a stellar core as in previous papers (e.g. Duffin fc Pudritzl2009t Baneriee et al 



20041). by e mbedding a slightly over-critical 1 M© Bonnor-Ebert sphere ( Bonnoi 1956 



Eberti 119551 ) in a low density environment. We add to this density distribution a 10% 
ove r-density to ensure collapse and a 10% m = 2 perturbation to break symmetry (see 
e.g. Duffin fc Pudritj|2009t LBaneriee et al. 2004 ). The background is an isothermal, low 
density environment in pressure equilibrium with the sphere (the density is set by choos- 
ing a background that is 10 times warmer than the sphere). The box is roughly 10 
times the size of the BE radius (0.81 pc in these models). We add to the sphere uni- 
form rotation such that the ratio of ro tational and gravitational ene rgy is moderate 
(Aot=0.046, similar to the value used in Hennebelle fc Fromane ( 20081 )) and a constant 
/? p iasma = 2ci? / 1 v\ = 46.01. In this model, one can relate the mass to flux ratio /i/^o, 
where fio is the critical mass to flux of /iq = , by fi/fio — 0.74c s /ua = 3.5. 

We use sink radii of 12.7 AU for longer runs, and radii of 3.2 AU to test the effect of 
sink particle size on the result. The accretion radius of a sink corresponds to 2.5 cells 
at the highest refinement level, and the critical gas density (beyond which gas can be 
accreted into particles) corresponds to the Jeans' density at the core temperature (20 K) 
of these cells. This gives p acc — 3.69 x 10~ 12 g cm~ 3 and p acc — 5.91 x 10~ u g cm~ 3 
for 12.7 and 3.2 AU sinks respectively. Each Jeans' length is refined by at least 8 cells 
and de-refined for more than 32 cells. Our customized version of the FLASH AMR code 
is described in previous wor k (|Frvxell et al. 1 120001 iBaneriee et al.ll2006l iDuffin fc Pudritj 
2008t iFederrath et al.|[2010h . 



3. The Early Collapse Phase 

The advantage of not using a sink particle is that we properly resolve the gas in the 
collapse. The disadvantage is that we are limited to pre-Class times (« 10 5 yr). We 
compare the following early collapse models: i) ambipolar diffusion, ii) ideal MHD and 
iii) hydrodynamics (e.g. B — 0). The rotational properties of these early structures are 
shown in Figure la. We compare surface density averaged values of v^/vr and w^/^Kepier 
to measure the degree of rotation in these early collapsed structures. The comparison is 
done at the maximum common central surface density (S 2 ). This is limited primarily by 
the ambipolar diffusion model which has the most constrained timestep. More evolved 
versions of ideal MHD and hydro models are included as thin lines in the graph, hinting 
at the evolution of the possible structures. The inward velocity has not halted (e.g. the 
flow is sub-Keplerian) in these early structures, however we are indeed seeing flattened 
rotationally dominated accretion discs (with associated outflows) being generated with 
sizes of 6-7 AU, facilitated by the suppression of gravitational instabilities. Bars produced 
in the hydrodynamic model have prevented a similar structure from forming by efficiently 
redistributing angular momentum. 

To demonstrate the suppression of fragmentation further, we ran the simulation with 
10 times more rotational energy in order to study early fragmentation (/3 ro t = 0.74, 
corresponding to high v alues seen in simulations of core formation in a turbulent medium 
(|Tillev fc Pudritzll2007l) ). ] ndeed our hydrodynamic collapse produces a wide binary with 
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Figur e 1. The early to late evolution of the protostellar core with ideal MHD. In a), (adapted 
from lDuffin fc Pudritd[2009l ), the rotation properties v^/vr and v^/vKepier of the early collapse 
(no sink particle). In b), contours of density (purple or medium grey, 3.3 x 10~ 17 g cm -3 and 
black, 1.33 x 10~ 15 g cm -3 at 2000 and 100 AU radii respectively) and v z (yellow or light grey, 
lkm/s) at the end of the simulation (with 3.2 AU sinks). In c), rotational properties of disk 
at the end of the simulation (with 3.2 AU sink). In d), mass evolution of different sink particle 
models, from left to right: 12 AU sink particle and a polytropic equation of state, 12 AU sink 
particle and 3.2 AU sink particle both with molecular cooling. 



a separation of about 1000 AU. The ideal MHD model suppresses the fragmentation , 
as documented in previous studies ( Price fc Bate 2007 : iHennebelle fc Tevssierl [20081 ) . 
resulting in a large bar and a central condensed structure. Meanwhile, the ambipolar 
diffusion model produces an intermediate result, a bar that has fragmented to form a 
binary. It perhaps through ambipolar diffusio n that any "fragmenta tion crisis" can be 
solved. These results are discussed further in ( Puffin fc Pudritdl2009l) . 



4. The Later Evolution: Sink Particles 

Using sink particles and ideal MHD, we were able to take the next step and evolve the 
collapse to much later times. The result is shown in Figure lb, wherein evolution over an 
additional 30 kyr is achieved. The sink particle mass evolution is shown in Figure Id for 
12 AU sink particles (with m olecular cooling or polytropic equation of state similar to 
that of Machida et al.l (|2010l) ) and a smaller 3.2 AU sink particle with molecular cooling 
( Baneriee et al.l 20061 ). By the end of the smaller sink collapse, the sink particle is O.2M 
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and the outflow has grown to 10 4 AU above the mid-plane of the disk. We are able to 
run the simulations of larger sinks further, and indeed these particles end up with nearly 
all of the core mass (80-90%). This is much higher than estimates of theoretical models 
seeki ng to relate the Core M ass Function (CMF) to the Initial Mass Function (IMF) 
(e.g. iMatzner fc McKed |2000|) . These results suggest that not much mass has in fact 
been cleared out by the outflow. This may be due to the fact that most mass must settle 
into the accretion disk before the outflow is launched, leaving less mass available to be 
cleared. 

The accretion disk is represented by the ratio of v^,/v r in Figure lc, showing a 2000 
AU accretion disk, which appears to be dissipating by the end of the simulation (as seen 
in Figure Id). There are two types of outflows, one lower speed (< 1 km/s) and a cen- 
tral, higher speed centrifugally driven wind. The inner disk and central component of the 
outflow are warped and precessing, as shown in Figure lb. The appearance of precessing 
warped discs and their outflow is extreme ly interes ting and is related to the back reac- 
tion of the MHD outflow on the disk fe.g. lLaill2003l) . Many of the qualitative properties 
of molecular outflows, including j et precession, dumpiness and an onion layered veloc- 
ity structure (jPudritz et al.1 120071 ) are seen in our simulation, and occur naturally as a 



consequence of solving the gravito-magnetohydrodynamic equations. 



5. Summary 

In the early stages, accretion discs are small (< 10 AU), massive, flattened, rotation- 
ally dominated, and are held together by the magnetic field. As the collapse continues 
over an additional 10 5 yr, the accretion disc and its associated outflow grow in size. The 
outflow torques and warps the disk leading to disc and jet outflow and precession respec- 
tively. Most material will be accreted by the star raising issues concerning the extent of 
protostellar feedback on stellar masses. 
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